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Radio -Wave Scattering by Tropospheric Irregularities 

Albert D. Wheelon ' 

(April L6, L959) 

The subject of radio-wave scattering by turbulent irregularities in the tropospheric 
index of refraction is reviewed. Descriptions of the turbulent medium are considered first, 
together with a status report on physical theories for the spectrum of irregularities. Phase 
and amplitude scintillations induced on electromagnetic waves propagated along Line-of- 
sight paths are discussed next. The complementary problems of signal statistics and elec- 
tromagnetic-propagation calculations are summarized and compared with available data. 
Both the geometrical optics and wave theory approaches are discussed. The theory of 
propagation beyond the horizon by scattering from the same irregularities is then described. 
A current review of the predictions for received power, signal fluctuations, and antenna 
effects is given. The paper is primarily a review of the essential and auxiliary predictions 
of scatter theory, but also contains a considerable amount of unpublished research by the 
author. 



1. Introduction 

The existence of random irregularities in the atmos- 
pheric index of refraction has been well established 
by direct measurements with microwave refractom- 
eters [1, 2]. 2 The influence of these refractive 
irregularities on electromagnetic waves is confirmed 
by the random fading of microwave signals received 
beyond the optical horizon. Stochastic phase and 
amplitude variations imposed on electromagnetic 
waves propagated along fine-of-sight paths confirm 
this cause-effect relationship. This is not to deny 
the influence of other mechanisms (layers, ducts, 
et al.) which may also contribute to scatter-field 
strengths. It does certify that the electromagnetic- 
ally predicted relation between (turbulent) refractive 
variations and signal characteristics is indeed es- 
tablished. 

The entire subject of electromagnetic propagation 
through turbulent refractive irregularities has re- 
ceived considerable theoretical and experimental 
attention during the post-war era. This attention 
is due to the enormous importance of "scatter 
propagation" communication systems on the one 
hand, and a basic limitation imposed on radio track- 
ing and guidance systems for space vehicles on the 
other. The subject also has considerable scientific 
interest in its own right, particularly insofar as it 
unites electromagnetic theory, turbulence theory, 
and probability in an essential mixture. Once es- 
tablished, the influence of atmospheric irregularities 
on electromagnetic waves presents a valuable tool 
for studying turbulence at very large Reynolds 
numbers by simple radio means. 

It now seems appropriate to try to establish how 
far the considerable volume of original research 
has brought the subject. This is particularly im- 
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portant since the basic concept is periodically brought 
under attacks of various sorts, and it is desirable 
to try to set forth the essential and auxiliary pre- 
dictions of the basic theory in a coherent, unam- 
biguous fashion. More important, however, is the 
need to present a balanced status report on the 
theoretical predictions, against which accumulating 
experimental results can be ranged to establish the 
true role of this mechanism in an objective manner. 

The present paper is primarily a review of the- 
oretical research on line-of-sight and scatter prop- 
agation in the United States, with a short status 
report on turbulence theories of atmospheric irregu- 
larities. It is not a genuine review, in that it leans 
rather heavily on the author's own research as a 
basis for unifying diverse results. Furthermore 4 , a 
substantial amount of unpublished material is 
included herein. The report goes one step further 
and points out a number of outstanding problems 
which are solved in principle, but not in detail. A 
detailed, balanced comparison of these theoretical 
results with experimental data was not attempted 
for this broad field, although occasional references to 
significant results and prior comparisons are indicated 
by way of stimulating this crucial testing. 

The substance of this report can be summarized 
briefly. In so far as the single scattering (Born) 
approximation and/or geometrical optics are valid 
representations of the electromagnetic response of 
propagating waves to refractive irregularities, the 
entire subject is solved in principle. The geo- 
metrical-propagation effects can be computed in 
any desired detail (or at least be reduced to integrals), 
independent of the turbulence model employed. 
This important separation is achieved by using the 
spectral representation for tin 4 space correlation of 
refractive irregularities, which allows any spectrum 
model to be combined with any propagation model. 
The more difficult electromagnetic problem of 
multiple scattering has not been considered ade- 
quately and is not discussed here. 
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2. Description of the Turbulent Medium 
2.1. Statistical Characterization 

The tropospheric index of refraction varies from 
point-to-point in a random manner. To describe 
this effect, one decomposes the dielectric constant 
(e=n 2 ) into its mean value and a small stochastic 
component. 



e(r, t)=e +Ae(r, t) . 



(2.1.1) 



It is this quantity which enters into the basic elec- 
tromagnetic equations. € is usually a function of 
altitude (h) only, and its gradient appears in certain 

physical theories of turbulent irregularities. The 

-> 
stochastic function Ae(r,0 is treated as a stationary, 
zero-mean random process in space and time, and 
is usually assumed to be spatially homogeneous and 
isotropic. If brackets denote ensemble, space, or 
time averages (which are equivalent for a statistically 
stationary, homogeneous atmosphere) the first two 
moments of Ae are: 



<Ae(r, t))=0 



(2.1. 



(Ac(r, t)Ae(r+R, t)) = (Ae 2 )C(\R\). (2.1.3) 

The intensity of dielectric fluctuations (Ae 2 ) is 
measured experimentally [1, 2] as several parts in 
10~ 12 . The space-correlation function C(R) falls 
to 1/g in a distance of several hundred feet, which is 
related to the scale length l Q of the irregularities. Little 
attention has been paid to higher moments of Ae, 
except to note that they are completely specified in 
terms of C(R) if the Ae form a Gaussian random 
process. 



2.2. Space Correlation Models 

The space-correlation function C(R) is the funda- 
mental characterization of the turbulence. First - 
order propagation quantities depend upon integrals 
of C(R), and one must assume an analytic model to 
perform the relevant calculations. The correlation 
must have the form 



C(R)=f(R/l ), 



(2.2.1) 



where f(x) is the Fourier transform of a non-negative 
function. l is called the scale length or average 
size of the blobs in the hierarchy of decaying eddies. 
The popular models iorJ(x) are plotted together in 
figure 1. 
1 . Bessel: j(x) = xK x (x) 

This correlation model was devised to predict: 
(a) The observed linear-wavelength dependence of 
average power for tropospheric scatter propagation; 
and (b) certain features of line-of-sight propagation 
experiments. It has zero slope at the origin, and 
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Figure 1. — Plot of four common space correlation models for 
isotropic atmospheric irregularities. 



is suggested by the theory of turbulent mixing-in- 
gradient. 

2. Exponential: f(x)=e~ x 

This function has been widely used to study scatter 
propagation. However, its cusplike behavicr (finite 
slope) at the origin implies infinite values for angle-of- 
arrival errors and other phase-derivative processes. 
The exponential model is possibly suggested by the 
physical theory of turbulent pressure fluctuations. 

3 . Gaussian : j(x) = e~ x 

This is widely used for line-of-sight calculations 
because it can be integrated easily. It does not 
predict the results of tropospheric scatter experiments 
adequately and apparently has no foundation in 
physical theories of turbulence. 

4. Cauchy: f(x) = [l+x 2 ]~ 2 

This model has been used for line-of-sight ray 
theory calculations, but is not founded in physical 
theories of turbulent irregularities. 

It should be emphasized that any correlation function 
is only a model of the atmospheric irregularities, and 
must be carefully evaluated in terms of its ability to 
predict the results of experiments. 

It is possible to generalize the foregoing to include 
anisotropic and inhomogeneous turbulent structures. 
Anisotropic irregularities can be described by intro- 
ducing different scale lengths in orthogonal direc- 
tions, viz, 



0(B)- 



f (\M + S") (2 - 22) 



and such effects are apparently relevant in the 
troposphere. 

Refractometer experiments show that the intensity 
of dielectric fluctuations decreases with altitude [1]. 
One can describe an inhomogeneous field of this sort 
with a height-dependent (Ae 2 ) in eq (2.1.3), so long 
as the intensity does not change significantly in the 
several scale lengths l Q required for the (space) 
averages. 
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A more difficult problem is encountered with the 
assumption of stationarity. There now seems to be 
experimental evidence [3] that the dielectric varia- 
tions either: (1) do not form a stationary process; or 
(2) that there is a great deal of power in the low- 
frequency part of the spectrum of Ae, which appears 
as trend when viewed over a short span of time or 
space. Bo tli explanations are probably partially 
correct, with a rapidly fading stationary process 
superimposed on long-period (diurnal, seasonal) 
harmonic components. Unfortunately, there is no 
adequate theory of nonstationary processes available 
to describe such effects. 

2.3. Turbulence Spectrum Method 

An equivalent characterization of a turbulent 

medium is given by the spectrum of irregularities 

-» 
S(k) [4, 5], which is the Fourier decomposition of 
the space correlation function, 



(^)C(K) 



~8tt 3 J 



d 3 ke- ikR S(k). (2.3.1) 



This representation has the advantge of exhibiting 
an explicit R-dependence for the correlation without 
committing the study to a particular turbulence model. 

Propagation calculations can be performed on the 

-> -» 
generic term exp [ — ik-Ii] by reversing the orders 
of integration. The physics of the atmospheric 

irregularities re-enters the problem eventually in the 

-> 
spectrum S(k), which must be introduced to com- 
plete the wavenumber (k) integration. If the 

-» -> 

medium is isotropic, S(k) depends only on |/r|, and 
the space correlation becomes: 

(A^C(h>)=^j^ dkkS(k) 5**£ (2.3.2) 



The mean square fluctuation at a point, 
** Jo 



(2.3.3) 



indicates that the spectrum represents the ability of 
each blob-size group (wavenumber range) to produce 
irregularities. 

Equation (2.3.1) may be inverted to establish the 
spectrum which corresponds to each of the correla- 
tion functions discussed earlier. Table 1 indicates 

Table 1. Correspondence of turbulence spectrum and space cor- 
relation functions for several models of atmospheric irregularities 



correlation 

C (R) 

Turbulence 

spectrum 

8 (k) 



Bessel 






«>**' 



[l+fc2$5/2 



Exponential 



e -R/l Q 



8t 



U+* 2 ® 2 



Gaussian 



- R2 ill 



x3/2(A C 2)/3 e 



Cauchy 



[■+?]" 



i-2(Ae2)^ -/c/ 



the spectra for the Bessel, Exponential, Gaussian 
and Cauchy models. The small It behavior of the 
correlation is determined by the large k range of the 
spectrum and vice versa. 

2.4. Physical Theories of Dielectric Irregularities 

The index of refraction in the (nonionized) tropo- 
sphere and stratosphere is independent of radio- 
frequency from 10 to 10, 000 Mc. Variations of the 
refractive index from unity are measured as several 
hundred parts per million, and are related to standard 
meteorological parameters by the empirical formula: 



n-} = 



1 77.0 J+3.73 10 5 ~ J-10- 6 (2.4.1) 



T = air temperature is degrees Kelvin (absolute), 

P = total air pressure in millibars, and 

e = partial water-vapor pressure in millibars. 

Dielectric variations (e=n 2 ) are thus related to tem- 
perature, pressure 1 and water-vapor fluctuations by 



Ae = 2 10- 6 < 



-Ar(77.6^+7.4(>10 5 ^ 

, A d P 7A] \ . A A73 10 ^ 
+ APl cr ) + Ae(- T ^) 



* (2.4.2) 



Dielectric variations in the lower troposphere are 
governed principally by water-vapor irregularities, 
and one is dealing with the irregularities established 
in a passive scalar by turbulent movements of the 
random atmospheric velocity field. This process is 
characterized by the turbulent mixing theories de- 
scribed below. 

In the stratosphere there is a negligible amount 
of water vapor, and dielectric variations are produced 
almost exclusively by thermodynamic fluctuations: 
AT 7 and AP. The primary mechanism for producing 
irregularities in the stratosphere is also probably 
turbulent mixing by rotating eddies. However, one 
must now recognize thermodynamic expansion (per- 
haps adiabatic [6]) of the gas as the eddies transfer 
parcels of air between different altitude 1 levels, cor- 
responding to different ambient temperatures and 
pressures. This expansion correction is not domi- 
nant in the lower troposphere, but probably plays 
a central role in the propagation of scatter signals 
to distances greater than 700 km by common volume 
elements in the stratosphere [6]. 

The spectrum method makes direct contact with 
physical theories or turbulent irregularities which 
predict S(k) in certain wavenumber ranges. These 
theories identify three distinct ranges of the spec- 
trum, as shown in figure 2. 

1. 0<£<& : The input or blob-creating range 
corresponds to very large scale lengths, which are 
probably anisotropic. The upper end, k , of this 
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Figure 2. — Typical spectrum of dielectric irregularities plotted 
versus wavenumber, showing various turbulence ranges. 



range is approximately the reciprocal of the correla- 
tion or scale length / discussed earlier, which is of 
the order of 100 m in the troposphere. Since we do 
not know (in detail) how turbulence is created from 
a laminar flow, turbulence theories carefully avoid 
this range. It is important to note that the spectra 
listed in table 1 do make definite statements about 
this range, so the large R behavior of the correspond- 
ing correlation models must be regarded as extra- 
polation or pure guesswork. 

2. k d <^k<Ck s : The inertial subrange is character- 
ized by redistribution of the turbulent "energy" 
toward higher wavenumbers and represents the 
progressive subdivision of eddies in the turbulent 
structure. Dimensional analysis can be used here 
with considerable confidence. Tropospheric-scatter 
propagation depends principally on this inertial 
range. 

3. ks<^k<^ oo : The dissipation range is character- 
ized by a sharp drop in turbulent activity, due to the 
destructive action of viscosity and diffusion. In the 
troposphere, l s =k~ l is of the order of millimeters. 
We do not consider the dissipation range, since 
aperture and recorder smoothing effectively elimi- 
nate the contributions of such tiny blobs in the 
troposphere. 3 

Three turbulence theories have been developed 
which predict S(k) in the inertial range: 

a. Pressure Fluctuations 

Fluctuations of the velocity field induce pressure 
variations according to Bernoulli's law. If variation 
of pressure, density, and water-vapor content are 
related linearly, the spectrum of dielectric fluctua- 
tions is given by [7], 

S(k)=QJ §^, (k>h) (2.4.3) 

where v Q is the speed of the largest blobs and c the 

3 The dissipation is apparently quite important for ionospheric-scatter propa- 
gation at vhf. 



local speed of sound. Most writers now agree that 
this mechanism is unimportant for the lower tropo- 
sphere [8]. 

b. Obukhov's Mixing Theory 

In this theory, an external source is imagined to 
feed fluctuations into the spectrum at the largest size 
(&o). Subsequent redistribution down the inertial 
range is attributed to turbulent convection, for which 
dimensional arguments give [9, 10, 11] 



S(k) = (AS) 



ftp 



(k>h) 



(2.4.4) 



If the external source is identified only with the 
mixing of a gradient by the largest blobs, the in- 
tensity is given by 



(Ae 2 )=fc - 



\dh)' 



(2.4.5) 



which is supplied at the input wavenumber k Q . 
Adiabatic expansions of the mixing eddies modify 
the ambient gradient in the stratosphere [6]. 

c. Mixing-in-Gradient Theory 

This approach considers an initial refractive-index 
gradient and how it is changed by turbulent con- 
vection at all wavenumbers. The rotating turbu- 
lent eddies transfer parcels of water vapor from low 
to high points on the profile and vice versa, and 
establish irregular contrasts thereby. The mech- 
anism for turbulent fluctuations appears explicitly in 
this theory at all wavenumbers and no external 
source is required. The inertial range is again 
amenable to unambiguous dimensional arguments 
[8, 12]. 

S(k)= (^Jw {k>ko) (2 - 4 - 6) 

and the result is independent of turbulent param- 
eters, except for restrictions on its range of validity. 
Adiabatic expansions of the eddies are again im- 
portant in the stratosphere [6]. In the range 
k^>k , we see that this model agrees with the spec- 
trum computed from the Bessel space-correlation 
model indicated in table 1. 

Both of the mixing theories predict that the 
intensity of turbulent fluctuations should be pro- 
portional to the local gradient of the inhomogeneous 
mean profile. A good correlation between the 
strength of scatter signals and the monthly mean 
refractive gradient has been observed experimentally 
on many paths [13]. Since strong scattered signals 
are associated with large dielectric fluctuations, 
this must be considered an important success for 
the gradient -dependent mixing mechanisms. On 
the other hand, observed absences of dielectric 
irregularities at inversion layers and other stable 
meteorological interfaces suggests that the mixing- 
in-gradient models may not be applicable to very 
sharp gradients. 
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2.5. Space-Time Correlations of Irregularities 

The averaged product of two dielectric fluctua- 
tions taken It different positions and times is of 
considerable importance in studying the time vari- 
ability of sigiiaj. amplitude and phase for both line- 
of-sight and beyond-the-horizon propagation. 



(Ae(r,t)^e(r+R,t + T))=(Ae*)C(R,T) (2.5.1) 

Temporal variations of Ae are associated with motion 
of the carrier medium, and two distinct types of 
atmospheric moticm are important : 

q. Drifting Convection 

-> 

A prevailing wind with constant velocity U 
would bear a frozen turbulent structure along intact 

and would translate the entire structure a distance 
-> 

U r in a time r. Since the atmosphere is assumed 
to be homogeneous, this would be just the same as 

measuring the space correlation between irregu- 

-> 
larities separated by a distance of U r at the same 
time. 

b. Random Motion 

The random motion associated with the con- 
tinuous eddy-breakup process rearranges the tur- 
bulent irregularities in time. This self-motion effect 
is probably isotropic for small (k) scale irregularities. 
The random motion would be apparent even if there 
were no drift velocity, or if one rode along on the 
prevailing winds. 

The spectrum method is particularly well-suited 
to describing the combined effecl of these two 
motions, in that their etl'eets can be explicitly 
separated in the integrand of the wavenumber 
integral representation of the 4 space-time correla- 
tion [4, 14]. 

{Ae(r,t)Ae(r+}i,t + r))^, ( d 3 kS(k)7 ik7 ^ Ur) rj(k 7 r) 

(2.5.2) 

For an isotropic spectrum of irregularities, this 
becomes : 



{Ae(r,t)Ae(r+R,t + r)) 



= * rdkkSifO Sm glf^ ,(*,r). (2.5.3) 
2T Jo \B+Ur\ 

The function rj(k,r) describes the time autocorrela- 
tion of fluctuations contained in a fixed wavenumber 
interval k, &nd is unity for zero time displacement r. 
Dimensional turbulence arguments indicate that 
r}(k,r) should depend only on the product 4 F &o /3 & 2/3 t 
in the inertial subrange [15, 16], although no func- 



tional dependence has yet been established. How- 
ever, detailed studies of the turbulent velocity 
field-time correlations over the entire spectrum 
indicate that the coefficient of t varies from F /3 to 
k[17]. 

Convective motion is probably the dominant 
effect for line-of-sight propagation, since eddy 
speeds F are usually much less than the prevailing 
wind speed U. 5 However, scatter propagation 
geometry discriminates against a horizontal drift 
velocity at the path midpoint and Do p pier shifts 
introduced by the random motion can dominate. 

3. Line-of -Sight Propagation 

3.1. Introduction 

Phase and amplitude instability on line-of-sight 
propagation paths is important in practice because 
of the limitations it imposes on interferometric track- 
ing schemes and frequency synchronization opera- 
tions. Before Launching into detailed discussions of 
problems, it is well to examine the two principal geo- 
metrical situations shown in figure 3. These prob- 
lems are idealized by omission of ground reflections, 
which are ignored since they do not contribute scin- 
tillations unless moving foliage or water surfaces are 
a major factor. 




&' 



P 



xy 



Primary Field 

Scattered Field Ray Path 

Figure 3. — Typical ray paths for idealized line-of-sight prob- 
lems: (a) Relay link immersed in an infinite turbulent me- 
dium, and (b) radio star viewed vertically through a turbulent 
mantel. 

a. Relay Link Problem 

Receiver and transmitter are both immersed in an 
infinite turbulent region a distance L apart. The 
transmitter emits spherical waves and the receiver is 
also assumed to be omnidirectional. 

b. Radio Star Problem 

Plane waves are assumed to fall on a semi-infinite 
turbulent medium with the receiver located a dis- 
tance L below the upper boundary. The turbulence 
may be graduated toward its upper boundary by 



* Vo and fco are the speed and wavenumber at the input stage of the turbulent 
velocity spectrum. 

5 Since the turbulent velocity field itself is a (small) partition of the kinetic 
energy in the laminar (drift) field. 
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using an inhomogeneous turbulence intensity (Ae 2 (r)). 
It is significant that the receiver always stands in 
the (tropospheric) turbulent medium for line-of-sight 
propagation, and blobs which are just on top of the 
receiver play an important role. 

3.2. Signal Statistics and Vector Voltage Diagrams 

The typical ray-path trajectories shown in figure 3 
indicate that the field which falls on the receiver is 
composed of two terms : (a) The unperturbed primary 
field E propagated along the line-of-sight, and (b) 
the incoherent superposition of waves scattered to 
the receiver by all volume elements in the scattering 
region V. The primary field induces a constant 
voltage A in the receiver, which serves as a con- 
venient phase reference. Because of the randomly 
varying multipath combinations, the total scattered 
wave induces statistical signal components x(t) and 
y(t) which are in and out-of-phase with respect to 
the constant vector A , as shown in the vector volt- 
age diagrams of figure 4. Experiments measure the 
composite amplitude A and the relative phase excur- 
sions a about the phase reference A . A typical 
phase record is shown in figure 4, together with the 
vector voltage diagrams appropriate to two displaced 
times ti and t 2 . 





a(t) 




Figure 4. — Typical phase record of randomly-varying line-of- 
sight signal showing the corresponding time-displaced vector 
voltage diagrams for a constant signal phis noise. 



The basic statistical problem is to estimate the 
probability distributions for amplitude A and phase 
a, and from these the relevant statistical estimates 
of the randomly varying quantities. Since the orthog- 
onal signal components x(t) and y(t) are induced by 
the scattering contributions from many independent 
blobs, they may be characterized by a Gaussian dis- 
tribution according to the Central Limit Theorem. 



P[x,y]dxdyt 



dxdy 



2lT(T x (Ty -yjl—p 

The moment parameters, 



exp— 



~x 2 .y 2 2pxy 



L 2(l-p 2 ) 



1 



(3.1.1) 



(x 2 ) = a 2 x , (y 2 ) = o- 2 v , (xy) = a x a yP 



should be determined from explicit propagation cal- 
culations. Previous statistical discussions have 
assumed, 

a^a^a 2 , and p = 0, (3.1.2) 

in order to simplify the analysis. We should like to 
point out, however, that such assumptions are not 
in agreement with available propagation calculations, 
which suggest rather that 



0"x+o"J = constant, p finite. 



(3.1.3) 



A new effort on the statistical problem is required for 
propagation work. 

For present purposes, we simply review available 
results derived on the assumptions (3.1.2). If the 
(A ) displaced cylindrical coordinates of figure 4 are 
inserted into eq (3.1.1), the joint amplitude-phase 
distribution becomes: 



t^/ a \ 7 a 7 AdAda 
F(A,a)dAda=— — $~ ex P 



-{ 



A 2 +A 2 -2AA c os a 

2(7 2 



} 



(3.1.4) 



The distribution of total voltage amplitudes A 
was established by Rice [18] from (3.1.4) by inte- 
grating the phase a from — w to t. 



5) 



P(^=^7 (^)exp-[^} (3.1 



From this expression one can calculate the ampli- 
tude moments. 



and 



+ 



(A 2 )=Al+2^ 



(0>m 



(3.1.6) 



(3.1.7) 



The probability that the amplitude exceeds a pre- 
scribed level a is obtained by integrating (3.1.5) 
from a to °° . Numerical tables of this probability 
have been prepared [19] and are reproduced on 
Rayleigh graph paper in figure 5. The addition of a 
constant vector evidently increases the average signal 
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Fir, tike 5. — Probability that the total amplitude exceeds a given 
level, measured relative to the amplitude of the con slant vector 
A , from [19]. 



level distribution, but reduces (lit 1 range of amplitude 
fading relative to A . 

The distribution of phase is derived from eq 
(3.1.4) bv integrating over all values of the ampli- 
tude A [20]. 

P(a)da=^ e~ A2 ° f2 ° 2 { 1 + J^ COS a e * ™ 2 ^ 



r j A Q cosa \"i1 



(3.1.8) 



Tliis distribution is plotted in figure 6 and indicates 
that the phase is evenly distributed when the con- 
stant vector Aq is very small compared with the 
rms multipath-scatter vector amplitude ^2a. When 
the constant vector is relatively large, Bremmer [20] 
proves that the phase is distributed about zero in a 
Gaussian manner, 



lim P(a)- 

<r«Ao 



A 



<jtJ2tt 



exp- 



{§:>' 



(3.1.9) 



Norton et al. [21] have also calculated the probability 
that the phase exceeds a given angle for various 
values of A /\ 2c The first and second moments of 
phase deviations have been evaluated numerically 
[22] using the result (3.1.8), and are reproduced in 
figure 7 as functions of 2a 2 / Al. 
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RATIO OF NOISE TO STEADY SIGNAL POWFR 



Figure 7. — Average phase variation and root mean square 
phase variations about constant references voltage A as func- 
tions of the amplitude ratio of signal-to-noise, from [22\. 



Bremmer [20] has considered the probability dis- 
tribution for time-displaced amplitudes (A u A 2 ) and 
phases (a u a 2 ) corresponding to figure 4. His de- 
velopment is based on special assumptions (see sec. 
4.1.) for the time correlation between x(t) and y(t), 
but is quite interesting in that he establishes phase 
and phase-rate distributions for the same instant 
by manipulation of the probability density for time- 
displaced phases. 
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3.3. Electromagnetic Scattering by Irregularities 

The stochastic dielectric constant (2.1.1) enters 
electromagnetic theory through the relation between 
the electric vector and displacement vector. 

zU[e +Ae(r,*)]i. (3.2.1) 

When the magnetic field is eliminated from Max- 
well's equations, the electric vector is found to 
satisfy [23], 



-» 1 ft 2 -» 



&*•'■} 



(3.2.2) 



Time derivatives here operate principally on E, since 
the field oscillates much more rapidly than the 

medium. The vector components of E are mixed 
in the right-hand term and give rise to depolarization 
effects. The magnitude of this term is small for line- 
of -sight microwave propagation and is usually 
omitted. (See, however, references [24] and [25].) 
With a harmonic time dependence (fc / =co/c=27r/X), 
one finds 



[V 2 +k}(e +Ae)]E=Q, 



(3.2.3) 



which is the basic propagation equation. 

It is not possible to solve (3.2.3) exactly, since 

-> 
Ae(r,t) is an unknown (stochastic) function of posi- 
tion. The Born approximation is widely used to 
describe single scattering by such fluctuations. This 
is essentially a three-dimensional iterative solution 

of the wave equation, which gives the total field in 

-» -> 
terms of the unperturbed wave E (r) as 



E(R) = E (R)-k 2 f d 3 rG(R, r)Ae(r, t)E (r) sinx 

(3.2.4) 

where V is the volume of irregularities illuminated 

by both the transmitter and receiver, X is the angle 

-> 
between E and the Poynting vector of the reradiated 

wave (see fig. 21) and is nearly 90 degrees for line-of- 

-» -^ 
sight propagation. G(R, r) is the free-space spherical 

Green's function connecting the scattering point 

-> -> 

r and receiver R. 



47r|#-r| 



(3.2.5) 



The first term in (3.2.4) represents the trans- 
mitted wave and the integral expression is the scat- 
tered wave. The total field induces voltage com- 
ponents in the receiver, as shown in figure 4. The 
primary field E sets the phase reference A . The 



integral term gives rise to the in- and out-of-phase 
random signal components x(t) and y(t), corre- 
sponding to the real and imaginary parts of the 
integral term in (3.2.4) respectively. The receiver 
measures the total voltage A and its phase a. It is 
inherent in the Born approximation that E s must 
be small compared with the mean signal E Q , so that 
instantaneous fluctuations of the total amplitude 
and phase are (consistently) approximated by: 

8A(t)^x(t)=kj T d 3 rAe(r, t) Re {G(R, r)B (r)}, 

(3.2.6) 



*(t)c 



y(t)_ kj 



A, 



K(R) 



i 



<PrAe(r,t) Im {G(R,r)E Q (r)}. 
(3.2.7) 



The variances of phase and amplitude fluctuations 
thus depend upon six-fold integrals of the space 
correlation function (2.1.3) 

It is common practice to use the geometrical 
optics approximation to describe line-of-sight propa- 
gation. This is equivalent to a one-dimensional 
(WKB) solution of the wave eq (3.2.3), which recog- 
nizes only the random speeding-up and slowing- 
down of the wave as it propagates along the ray 
path (s). The phase of the vector voltage diagram 
in figure 4 is given directly by this method, 



a= lT 



dsAe(s,t), 



(3.2.8) 



and is to be computed along the nominal (recti- 
linear) ray path s. This expression is always a good 
approximation to the volume integral (3.2.7), and 
it is generally used to study phase variations because 
of its simplicity. Geometrical optics also predicts 
amplitude variations of the received signal. These 
are due physically to bunching and diverging of the 
energy-bearing rays by random refractive bending ; 



8A 



=rJ>[S+$] <"*> 



where a is given by expression (3.2.8). This expres- 
sion is usually not a good approximation to the more 
precise wave theory result (3.2.6). 

The geometrical optics approximation is a special 
case of the single scattering wave theory of eq (3.2.4) 
[4, 5, 26]. To understand this correspondence, note 
that the scale length of the irregularities Z is ordin- 
arily much larger than the wavelength of the radia- 
tion. The blobs act like very shallow lenses and 
project almost all of the scattered radiation in a 
cone of angular opening X// about the line-of-sight. 
The effective scattering volume over a path length 
L, therefore, has a maximum width iX// as shown in 
figure 8. If this width is much less than the average 
blob size, the scattering volume consists of cylindrical 
plugs "cored" from successive blobs as in figure 8a. 
The propagation is then substantially one dimensional 
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Figure 8. Pictorial representation of the effective scatter volume 
for line-of -sight radio star problem: (a) Fresnel scattering, 
and (6) Fraunhofer scattering. 



and random variation of the phase speed of the wave 
is the measured effect. This corresponds to Fresnel 
scattering and is correctly treated by geometrical or 
ray optics. The opposite extreme occurs when 
LX/7r/ >>/ ; the receiver sees many whole blobs as 
in Cgure 8b. The problem is then necessarily three 
dimensional and requires the wave solution (3.2.4), 
corresponding to Fraunhofer scattering. 'Die scatter- 
ing parameter 



7=: 



Lk 

J 2 



(3.2.10) 



is evidently the important variable in these dis- 
criminations. One finds that the wave theory expres- 
sions (3.2.6) and (3.2.7) reduce exactly to the geomet- 
rical optical results (3.2.9) and (3.2.8) if y is less than 
unity. If 7 is greater than (or near) unity, how- 
ever, one must use the full-wave theory to compute 
amplitude scintillations. 

The single-scattering approximation (3.2.4) breaks 
down at large distances L and/or high frequencies, 
since the convergence of the Born series is related 
to the smallness of the mean square phase expres- 
sion [5], 



<a 2 >^7r 2 <Ae 2 > 



Ll 



(3.2.11) 



As the wave progresses deeper into the turbulence, 
scattering out of the main beam gradually dissipates 
the primary signal A in figure 4 and one is left only 
with the correspondingly increased random field 
E s . The method of Rytov [27] has been used to 
describe this "phase creep" zone, where the amplitude 
and phase fluctuations are not small. 



E(R)=E (R)exp- 



h 2 

_E (R) 



j: 



(PrG(I?,r)Ae(r,t)E (r) 



(3.2.12) 



This expression reduces to (3.2.4) in the limit of 
small signal variations, but it is not clear how this 
method accounts for multiple scatterings. It is 



significant that the integral term in the exponent of 
Rytov's expression is identical to the first Born term 
in eq (3.2.4). Difficult propagation calculations of 
the real and imaginary parts of this integral term 
computed for single scattering theory can thus be 
used in calculations based on (3.2.12), once the 
appropriate signal-statistical description is estab- 
lished. 

Feinstein [28] used the Kirchhoff-Huygens approxi- 
mation to study propagation of scalar waves through 
a succession of uncorrelated random slabs. His use 
of the second order phase approximation is equivalent 
to tracing the phase front along various ray paths, 
or neglecting diffraction over the slab thickness. 
The result can be expressed as a power series in 
<Ae 2 >, the first term of which is the geometrical- 
optics approximation. 

3.4. Propagation Calculations Based on Geometri- 
cal Optics 

The foregoing discussion presented general expres- 
sions for the phase and amplitude scintillations of 
waves which are propagated through random media. 
To compare theory and experiment, it is neces- 
sary to combine these propagation expressions with 
the descriptions of turbulent media given in part 
II. The results depend both on the propagation 
geometry and the model chosen to represent the 
irregularities. 

The calculation of phase variations is usually based 
on the geometrical optics approximation. The 
variance of phase fluctuations <Ca 2 ^> about the 
constant reference vector A can be computed by 
simply squaring and averaging the integral expres- 
sion (3.2.8). Using any of the correlation models 
described in section 2.2 gives the following simple 
result, 



2 > = g7r 2 <Ae 2 > 



Llo 



(3.3.1) 



where L is the line-of-sight path length and X the 
wavelength of the radiation employed, q is a con- 
stant of order unity which has been evaluated for 
each of the four correlation models. References to 
such detailed calculations are summarized in table 2, 
together with references to other phase and angle-of- 
arrival quantities. The reader should consult these 
references for specific results. 

The wavelength dependence of expression (3. 3.1) 
has been confirmed experimentally by simultaneous 
measurements [33, 34] of phase shifts on three fre- 
quencies, 100, 1,000, and 10,000 Mc, on line-of-sight 
paths to mountain tops. The path length (L) de- 
pendence is difficult to check experimentally, since 
paths of different lengths are apt to go through 
regions of the troposphere which do not have the 
same level of turbulent activity. Attempts to corre- 
late phase variations with airborne and ground-based 
refractometer measurements of <^Ae 2 > and l in the 
propagation paths have been moderately successful 
[3, 30]. 

Such results are conveniently summarized with the 
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Table 2. References for line-of-sight calculations of phase and 
amplitude fluctuations using geometrical optics for four space 
correlation models of the irregularities 



\^ Space 
\ correlation 
Calculated \ model 
quantity \ 


Bessel cKi(c) 


Exponential 


Gaussian 
e-/>2 


Cauchy 

[1+P2J-2 


Mean square phase 

Finite path corrections 
to<«2> 


30 
30 
31 


29,30 
30 
31 


29,30 
30 
31 


29,30 
30 


Aperature smoothing 
corrections to <C<x 2 ^> 




Finite data sample cor- 
rections to <a 2 !> 




Parallel path space cor- 
relation 




29 

29 

29,30 

29,30 
Diverges 

Diverges 

Diverges 

29 

29 

Diverges 


29 

29 

29,30 

29,30 
32 


29 


Rotated path space 
correlation 


30 

30 

30 
30 

Diverges 


29 


Time correlation of 
phase.. . --. 


29,30 


Power spectrum of 
phase 


29,30 


Phase fading rate 

Number of phase max- 
ima 


32 


Mean square angle-of- 
arrival 


29 
29 
29 


29 


Time correlation of an- 
gle-of-arrival 




29 


Spectrum of angle-of- 
arrival 




29 


Mean square ampli- 
tude variance 


Diverges 











spectrum method, suggested in section 2.3, where a 
model choice is delayed until after the propagation 
calculations have been performed. In terms of the 
blob-size spectrum S(k), the mean square phase 
fluctuation becomes [4] : 



<« 2 > 



2X 2 Jo 



dkkS(Jc), 



(3.3.2) 



which exhibits directly the dependence of the meas- 
ured (<a 2 >) on the controlled variables L and X. 
The spectrum method is particularly advantageous 
when computing corrections to this expression due to 
aperature smoothing, finite data sample lengths (T) 
and finite propagation path lengths (L). These 
effects can be analyzed exactly within the framework 
of geometrical optics, and produce appropriate fac- 
tors in the wavenumber integration which filter the 
spectrum S(k) preferentially. Analyzing the indi- 
vidual effects separately and multiplying them 
serially one finds [4]. 

x 



[|{ siTO _Hg*S}] 



X 



h^r^^ 21 ^} (3 - a - 3) 



This composite spectrum filtering is illustrated in 
figure 9. The first factor accounts for the masking 
of small blob (large wavenumber) contributions to 
<C<* 2 > due to smoothing by a receiving aerial of 



Filter Factors 




Figure 9. Typical spectrum showing preferential filtering 
factors: (1) finite path length, (2) finite data sample length, and 
(3) aperture smoothing corrections. 



radius a. The second factor describes finite path 
length correlations, and indicates that blobs larger 
than the total path length L cannot contribute 
effectively to signal fluctuations. The last factor 
describes the effect of a finite data sample length T 
in suppressing very low-frequency scintillations, 
assuming that they are due only to a convective 
speed U. This predicts that the rms phase variation 
should increase with sample length T, which is 
observed experimentally [3] ; although nonstationary 
trends could give the same result. References to 
explicit calculations of the individual effects for 
various spectral models (i.e., correlation functions) 
are given in table 2. 

The space correlation between phase variations 
measured on adjacent propagation paths has been 
calculated for a variety of turbulence models (see 
table 2) for both the rotated line-of-sight and 
parallel path configurations indicated in figure 10. 



•X l ■ r ~'~~. : '- 

\ ■ ' 

w — — ~ v& *& — 
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Figure 10. Geometry for space correlation of phase experiments: 
(a) Rotated paths with common transmitter and (b) parallel 
rays from source of infinity. 
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The spectrum method summarizes these results in 
succinct form [4]. 
Parallel rays : 

<a(x)a(x+D)>=^ 2 i "dkkSWJoikD), (3.3.4) 

-A JO 

Rotated rays: 

<«(*)«(*+0)>=|^J o dkkS(k)~j Q duJo(u). 

(3.3.5) 

These filter factors essentially eliminate the con- 
tribution of small blobs (large wavenumbers), which 
are not large enough to affect both paths simul- 
taneously. 

It is found that the Bessel model gives good 
agreement with spaced receiver experiments [33,30] 
for small baseline distances. At very large separa- 
tions, however, the experimental correlation falls 
more slowly than the theory predicts. This is 
tentatively ascribed to lack of spatial homogeneity 
and/or temporal stationarity. Since the correlation 
at great distances is due to the joint influence of very 
large blobs on the two propagation paths, the phase 
correlation should depend on the large R behavior 
of C(R), which is not tied to physical theories, as 
noted in section 2.4. 

The random variation of single path phase records 
with time is due to the action of drifting convection 
and turbulent self motion. The space-time correla- 
tion of dielectric fluctuations (2.5.3) permits one 
to calculate the time correlation of phase records 
quite generally for an arbitrary spectrum. If the 
convective wind blows normal to the propagation 
path [4], 



<a(t)a(f+r)> 



2X 2 Jo 



(IkkSWJolkUrMk^]. (3.3.6) 



The effect of self-motion is usually discarded (i.e., 
17= 1) ? since the turbulent speed (V ) is much less 
than the drift speed U. 

The experimental data on phase scintillations is 
usually expressed in terms of frequency spectra of 
the phase records. For drift motion alone [5], 



W(f)- 



= 2X 2 



J** 



dkkS(k) 



[k 2 U 2 -(2wf) 2 ] 1/ 



(3.3.7) 



The mixing-in-gradient model (2.4.6) predicts that 
the frequency spectrum W(f) should vary as /~ 4 , 
while the Obukhov model (2.4.4) predicts an 
/~ 2 - 7 variation. The most recent data [3] indicates 
a frequency variation of /~ 2,8 from 1 cycle per hour 
to 10 cps for very low-level paths. This data argues 
strongly for the Obukhov Mixing Theory, unless 
self-motion effects play a dominant role. 

The average number of mean (zero) crossings and 
maxima can be expressed in terms of ratios of mo- 



ments of the frequency spectrum W(f), which in 
turn can be related to moments of the turbulent 
spectrum S(k) [4]. Filter factors at the high-wave- 
number end, such as that provided by aperture 
smoothing, are especially important if such results 
are to be compared with experimental data. 

Tracking devices which measure bearing angles 
ordinarily measure the normal to the instantaneous 
phase front by interferometric means. Angle-of- 
arrival errors are thus related to the space coherence 
(correlation) of phase scintillations along this phase 
front. The mean square angle of arrival is given by 
[4] 

< 5 02> = A r dWS(k), (3.3.8) 

l07rjo 

and has been evaluated explicitly in the references 
of table 2. The autocorrelation of 8d(t) and the 
corresponding frequency spectrum have also been 
discussed in general [4] and specific terms [29]. 

Amplitude fluctuations have not received as much 
attention as phase (see however, ref [36] and [37]) 
since 8A/A must be small compared with (a)rm 8 
for the geometrical optics expression (3.2.9) to 
be valid. In terms of the spectrum, the percentile 
amplitude fluctuation is given by 



8A 

A 



> p f 



dkk 5 S(k), 



(3.3.9) 



which depends critically on the behavior of the cor- 
relation model near the origin (i. e., k=co). Only 
the Gaussian and Cauchy correlation models pre- 
dict finite results, viz, 



bA 

A, 



''\=g(Ae 2 )^ 



(3.3.10) 



where g is a constant of order unity. This distance 
variation and absence of a frequency dependence 
does not agree with most experiments [33, 38], so 
one must seek a more satisfactory solution in the 
wave theories. 



3.5. Propagation Calculations Based on Wave 
Theory 

When the scattering process is necessarily three- 
dimensional in the sense of figure 8, one must use 
the Born expressions (3.2.6) and (3.2.7) for line-of- 
sight amplitude and phase. The variance of phase 
and amplitude depend upon six-fold integrals of the 
space correlation function (2.1.3). Such calcula- 
tions have been performed primarily with the 
Gaussian function, which has important analytical- 
simplification properties, although even then one 
must exploit approximations. 

The spectrum method has simplified the evalua- 
tion of these wave theory expressions enormously. 
Using the method described previously [4], one 
can establish the following general expressions for 
the radio-star problem : 
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where the fundamental scattering filter functions are 
defined by 



$(z) = 



¥(z) = 



2(l+z 2 ) 
a; 



2(l+x 2 ) 






(3.4.2) 



(3.4.3) 



The corresponding filter factors in the spectral inte- 
grals of (3.4.1) are plotted together in figure 11. 
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Figure 11. Fundamental scattering filter factors for wave 
theory treatments of line-of-sight propagation. 



These functions depend upon the scattering param- 
eter y=k 2 L\ for each blob-size group in the dis- 
tribution described by S(k) . Those blobs in S(k) 
which scatter in the Fraunhofer zone (t»1) 
produce equal phase and amplitude variations, each 
one-half as great as the ray theory result (3.3.2). 
The Fresnel regime (7-Cl) gives the geometrical 
optics phase result (3.3.2) directly. Using the 

small argument behavior $(ar)~~j the geometrical 

optics result for amplitude fluctuations (3.3.9) is 
regained. 

Explicit expressions for the phase and amplitude 
correlations can be established by introducing a 
particular spectrum model from table 1 into (3.4.1) 
and performing the single remaining wavenumber 
integration. Barrows [35] has evaluated these in- 
tegrals for the Exponential, Gaussian, and Bessel 
models, and finds 



hA 
A 



(n, 2 \- 


2 7 r 2 <Ae i 


>Ll 


\a ) 


X 2 




\y- 


2*- 2 <Ae< 


>Ll 


X 2 




SA\_ 


_2x 2 <Ae 


2 >Lh 



F 






[< 



<m 



Ao/ 



G (^Il) 



(3.4.4) 



where the scattering parameter-dependent functions 
have approximately the same form as the filtering 
functions &(x) and $(x). The results obtained in 
this way are substantially the same as those ob- 
tained by Fannin [42], Mintzer [6], Obukhov [39], 
and Chernov [40] by the more laborious direct 
calculation of six-fold volume integrals. Since their 
results differ in details of propagation geometry, 
space correlation models, and scattering theory, 
we have constructed table 3 to summarize the 
essential assumptions and results of each paper. 

The derivation of expressions (3.4.1) given in 
reference [4] does not lend itself to calculating space 
and time correlation effects. A new method was 
devised which combines the integral expressions 
(3.2.6) and (3.2.7) with the original representation 



Table 3. Explicit calculations of phase and amplitude fluctuation quantities using three dimensional wave theory 



Author 


Refer- 
ence 


Scattering theory 


Propagation model 


Space correlation 
model 


Phase 
variation 


Ampli- 
tude 
variance 


Phase- 
ampli- 
tude cor- 
relation 


Space 
correla- 
tion 


Time 
correla- 
tion 


Others 


Obukhov . _. 


39 
40 
41 

6 
42 
26 
23 
35 

5 


Rytov _ 


Radio star 

do 


Gauss 


7 
V 

V 

V 
V 
V 


1/ 


V 

J 


7 
/ 


V 
V 




Chernov. _. _. 


do 


do 




Tatarski___ _„ __ 




_. do 


Exp... 




Mintzer. _. __ 


Born. 


Relay link 

do 

Radio star 


Gauss. . . -. 


y 


Fannin. . 


do 

do 

do 


do 




Ellison.. ._ _ . 


____ do 




Wheelon 


.... do 


do 




Barrows,.- 


do 


do 


Bessel Gauss 

General. . 




Wheelon _.. _. 


do 


do 
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(2.3.1) of the space correlation function, rather than 
the specialized form of eq (2.3.2). For example, 



'J; 



(a 2 )=kj d?rlm[G{ll,r)E i 



itoif 



cPr'Im[G(R,r')E (.r')} 



X~J^kS(k) 



exp ih(r—r f ) 



=p, (<PkS(k) I ( <fa^ r Im[G(Rj)Eo(r)] 



(3.4.5) 



Barrows [35] has found that the volume integral 
can be evaluated rigorously for either the radio star 
or relay link problem. In principle one can include 
the influence of ground reflections indicated in figure 
12, by simply redefining the Green's function G(R,r) 
and incident field E (r) in eq (3.4.5) to satisfy the 
boundary conditions at the earth's surface. This 
can be accomplished analytically for a perfectly 
conducting earth by simply adding image sources 
and image scattering blobs to represent the ground- 
bounce paths. 
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Fir.uRE 12. Typical ray paths for line-ofsight propagation 
including ground reflected effects: (a) Relay link & (b) radio 
star problems. 

The new representation (3.4.5) is well adapted to 
describing time autocorrelations of the phase. Using 
the space-time correlation spectral representation 
(2.5.2), we find 



(a(t)a(t- 



r))=^jd*kS(kMk,T)e- 



X 



f cPre** 



Im [<?(i?,r)S (r)] 



(3.4.6) 



which uses precisely the same geometrical propaga- 
tion integrals required for the phase variance calcula- 
tion in (3.4.5). Space correlations follow directly 

-> 

if one replaces the drift displacement Ur with the 

-> 
receiver separation D, providing the appropriate dis- 
placed coordinate system indicated by figure 13 is 
employed. If the wavenumber integration is per- 
formed in polar coordinates, it is found that the 
geometrical integrals do not depend on the (k) 
azimuth angle, and therefore, 




Figure 13. Coordinates for calculation of phase correlation 
between adjacent receivers. 



(a(x)a(x+D))= 1 ^ 2 f dkk 2 S(k) ( 'dd sin 8j {kD sin 6) 
■">*" Jo Jo 



xl ( (Prelim [G(R, r)£i(rj] 



(3.4.7) 



The effect of aperture smoothing with a circular 
parabolic dish of radius a can be evaluated by the 
same integration technique used with the less ac- 
curate geometrical optics solution [4, 31]. 



I 2 \ 






IGtt 2 



II 



X (Pre*** Im[Gyi,r)Eo(r)]\ 



(3.4.8) 



Even the difficult question of polarization dependence 
can be discussed conveniently in this framework, if 
a tensor Green's function is used in the Born ap- 
proximation [24, 25]. In this way one can build up 
theoretical estimates of considerable interest to 
measurement programs. The method keeps the 
geometrical and turbulence problems separated until 
the last possible step, so as to allow each to develop 
as fully as possible. 

3.6. Physical Interpretation of Scattering 

It is valuable to pause briefly to re-examine the 
mathematical expressions (3.2.4) introduced by the 
Born approximation. Figure 14a indicates a typical 
singly-scattered ray trajectory for the radio-star 
problem. The scattered ray path concept may be 
identified with terms in the integrand of (3.2.4). 

The incident plane wave E (r) propagates recti- 

linearly to the scattering volume element d?r. The 

-» 
dielectric fluctuations Ae(r,t) scatters a small fraction 
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SINGLE SCATTERING 




DOUBLE SCATTERING 



Figure 14. Ray paths for single and double scattering of 
plane wave by a semi-infinite region. 



of the wave E (r) in all directions. (This reradiation 
is strongly concentrated in the forward direction, 
since the refractive variations act like very weak, 
overlapping lenses with radii of curvature l much 
larger than the radiation wave length X) . A portion 
of the scattered energy is propogated along the line 

connecting the receiver R and scatterer r. The 
propagation of this component is described by the 
Green's function (3.2.5), which characterizes the 
spherical expansion of the (unit) scattered wave. 
The total scattered field is simply the summation 
of all such multipath contributions, as given by the 
integration (3.2.4). These scattered waves are out 

of phase with respect to the primary signal E pro- 
pagated along the line-of-sight, in virtue of their longer 
path lengths. Random fluctuations of the strength 
of these scattered components (i.e., Ac) give rise to 
the characteristic fading suggested by figure 4. 

The typical ray path for double scattering pro- 
cesses indicated in figure 14b can be compared with 
the second Born approximation 



E(R) = E (R)-k 2 f f d z rG(R,r)Ae 



(r,t)E (r) 



+kj 



j>j". 



d z r'G(R,r)Ae(r,t)G(r,r') 



Ae(r' f t)Eo(r f ) (3.5.1) 

The double integral term has the following meaning: 

an incident wave E (r') falling on the first scattering 

-» 
element d?r' is scattered by Ae(r',£), and reradiated 

-*->-» -> 

to r by G(r,r'). The second element Ae(r,t) scatters 
the wave again and reradiates the result to the 



receiver via G(R, r) . The double integration simply 
sums the contributions over all possible pairs of 
scattering elements. 

This tracing of rays between scattering events is 
similar to Feynman's diagramatic technique for 
describing the space-time history of particles in 
Quantum Electrodynamics [43]. In treating the di- 
electric fluctuations Ae as a sequence of scattering 
events, however, one must give up the concept of a 
continuously refracting medium. Phase fluctuations 
are no longer regarded as the result of random 
changes in the velocity of propagation along the 
line-of-sight path. Rather, phase instability is looked 
upon as the random composition of many multipath 
components with the primary signal (see ref [44]). 

Our description has been something less than 

accurate thus far in that expansions like (3.5.1) 

-> 
imply more power at the observation point R than 
would have appeared in the absence of refractive 
fluctuations \E^\ 2 . We have evidently overlooked 
scattering out of the primary beam on the line-of- 
sight ray path to the receiver. A scattering event 
robs the beam of some energy at each point along 
this path. Such losses are just compensated by 
energy which is projected back into the receiver by 
off-axis blobs. Inspection of figure 15 shows that 
for each point to which the line point scatters energy 
out of the beam, there is a corresponding point on 
the wave normal which projects an equal amount of 
energy back into the receiver at the same angle. 
If the receiver is isotropic, the two effects just cancel. 




-o ;,f°"r s i \ G r in jx 




Figure 15. Singly-scattered ray paths for radio star problem 
indicating power balance between scattering losses out of the 
beam and projection of energy into the receiver by off-axis 
blobs. 



Energy can re-enter the primary ray only by a 
double scattering process, which changes its direction 
back to the original line-of-sight. This means that 
the power balance must be maintained at each level 
of iteration, and is related to the unitary relation- 
ships which must hold between the various terms in 
the Born series. It is important to observe that 
such conservation cannot be demanded of the fluctuat- 
ing field strength, but only of the averaged power. 
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4. Scatter Propagation 

The transmission of telephonic and even television 
signals well beyond the optical horizon at UHF fre- 
quencies is such an important component of world- 
wide communications today that it requires no 
historical introduction. Suffice it to say that weak 
but dependable communication signals are received 
out to distances of 1,000 km if sufficient transmitter 
power and receiver sensitivity are employed. The 
signal envelope is observed to fade in an irregular 
manner. This suggested that the propagation might 
be due to radiowave scattering by the refractive ir- 
regularities (turbulence) which were known to exist 
in the troposphere. The mathematical theory of 
this scattering was developed by Pekeris [45] and 
Booker and Gordon [46], using the single scattering 
approximation. 

The refractive scattering explanation of micro- 
wave propagation beyond the horizon is illustrated 
by figure 16. Turbulent blobs in the common volume 
of the receiver and transmitter are imagined to scat- 
ter minute fractions of the line-of-sight fields into 
the receiver. The case of omnidirectional aorials is 
shown in figure 16a, and indicates that the common 
volume is defined by the tangent planes to the trans- 
mitter and receiver. The incoherent addition of 
signal contributions throughout this volume gives 
rise to random orthogonal voltage components in the 
receiver, as suggested by the drawing. The narrow- 
beam case, shown in figure 16b, indicates that the 
corresponding scattering volume is controlled by the 
antenna patterns. 





Figure 16. Geometry for scatter propagation beyond the 
horizon: (a) Broad beam aerials, and (6) narrow beam aerials. 



A substantial amount of reliable experimental data 
was first published in 1955 [47]. The overall experi- 
mental situation at that time was ably summarized 
by Bullington's review paper [48]. Additional results 
have been published sporadically in succeeding years, 
and it is somewhat difficult to establish the precise 
state of affairs at this time. An up-to-date summary 
of theoretical research on the subject in the United 
States was presented to URSI in 1957 by Staras and 
the author [49]. The following section extends that 
report by describing a considerable amount of un- 
published work. We also attempt to correlate the 
theoretical predictions with experimental results 
where possible. 

4.1. Signal Statistics and Vector Voltage Diagram 

The electromagnetic field received beyond the 
horizon via the scatter mode induces a complex 
voltage 



x(t)+iy(t)=R(t) e** 



CO 



(4.1.1) 



in the receiver, as shown in figure 16. The orthog- 
onal signal components x and y result from the in 
and out-of-phase incoherent addition of waves 
scattered from many independent blobs in the com- 
mon volume. Application of the Central Limit 
Theorem to this scattering process by many blobs 
ensures that x(t) and y(t) are distributed in a Gaus- 
sian manner. The components change randomly 
with time because of the changing phase relationships 
of the scattering contributions from randomly mov- 
ing blobs. The processes x(t) and y(t) are probably 
stationary over brief observation periods, but do 
exhibit diurnal and seasonal variations which are 
quite predictable. 

The probability distribution for amplitude R(t) 
and phase <j>(t) at a fixed time t can be computed 
directly from the Gaussian distribution for x(t) and 
y(t) suggested above. It is usually assumed that 
both these components have equal variances a 2 in 
the far field scattering range of beyond-the-horizon 
propagation. 



P(x,y)dxdy= 



where 



dxdy 



2tt(7 2 Vi-P : 






2pxy *\ 



p 2 ) /' 



(x 2 )={y 2 )=° 2 



(xy} = 



Vp. 



Transformation to the polar coordinates of (4.1.1) 
yields the joint amplitude-phase distribution. 



P(R,4>)IMEd<l>= R( !f^- Z exp- 

JiTrcr y 1 — p 



f /Z 2 [l-psin 2 4> n 
\ 2<r 2 (l-p 2 ) / 



(4.1.2) 



Previous discussions of random signals have taken 
the cross-correlation (p) between the in and out-of- 
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phase components to be zero. This is justified for 
electrical shot noise which is a randomly phase- 
modulated signal. It is not justified for scatter 
signals, which are essentially amplitude-modulated. 
This cross correlation has been computed explicitly 
for line-of -sight propagation geometries [40], and is 
definitely not zero in that case. Although a similar 
calculation needs to be performed for scatter propa- 
gation geometries, it is quite possible that p has a 
finite value. 

The distribution of signal amplitudes is established 
by integrating (4.1.2) over all values of the phase, 

( — T tO 7r). 

(4.1.3) 



,If the cross correlation p were zero, this would reduce 
t o the Rayleigh distribution for R [50, 51]. The 
r oot mean square amplitude is independent of p, 



Rrms=^W) = ^2 



(4.1.4) 



and forms a convenient amplitude reference. 

The probability that the amplitude exceeds a given 
level r is obtained by integrating (4.1.3) from r to 
infinity, and is plotted in figure 17. The theoretical 
curves for different p are plotted on Rayleigh paper, 
so as to exhibit departures from the straight line plot 
of a simple Rayleigh distribution (p=0). This 
curve indicates that finite values for the cross 
correlation enhance large values of R relative to those 
of the Rayleigh distribution. The probability of 
observing very small signal values is correspondingly 
smaller than for the Rayleigh case. It is quite 
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PROBABILITY THAT SIGNAL EXCEEDS ORDINATE VALUE 

Figure 17. Probability that instantaneous signal amplitude 
exceeds a given level measured in terms of the RMS signal 
level for various cross correlations p. 



significant that this very departure has been ob- 
served experimentally on many scatter links [47, 52]. 
It should also be noted, however, that a signal level 
trend superimposed on the random amplitude process 
would produce the same general behavior for finite 
data samples. 

The possibility of a steady signal being admixed 
with the random scatter signal cannot be ignored, 
especially on short paths where the ground wave 
or reflecting layers may be important. The cor- 
responding statistical problem of a constant vector 
plus a random process was summarized in section 
3.1. It is apparent that the entire signal level distri- 
bution is raised by such an addition, which offers a 
third possible explanation for the non-Rayleigh 
character of the signal level distributions. 

The average signal level can be computed from. 
(4.1.3). 



<£>=< 



.vs 



r i[^Hi^)] ^ 



where E(x) is the complete elliptic function of the 
second kind. The correlation -dependent factor in 
square brackets varies by only ten percent over the 
entire range of p, so the Rayleigh form is an adequate 
expression. 

The distribution of signal phases can be estab- 
lished by averaging (4.1.2) over all possible ampli- 
tudes. 



p[m=P 



■y/l-f 



2tt 1— p sin 2<p 



(4.1.6) 



This is plotted in figure 18 for 0O<7r (since it is. 
symmetric about the origin) and various p, indi- 
cating the departure from a uniform phase distri- 
bution which is usually associated with a Rayleigh 
distribution of amplitudes. 

The foregoing description applies to the signal 
fluctuations observed during an interval ranging 
from several minutes to several hours. A sub- 
stantial amount of experimental data is presented 
in terms of distributions of hourly medians, which 
are found to be "log-normally" distributed [53]. 
McCrossen [54] used this experimental result to sug- 
gest a phenomenological theory of nonstationary sig- 
nal distributions for all time spans. The essential 
idea is to treat the signal variance in eq (4.1.3) as a 
long-term random variable. It is assumed that a is 
sufficiently constant over intervals required for the 
averaging of short-term fluctuations. The quantity 
which is observed to be log-normally distributed is 
the hourly median of received power, 

Z=P . 5 =(7 2 log 4 = <i? 2 > log 2, (4.1.7) 



where p is taken to be zero for 
distribution for Z is thus 



convenience. 



The 



P(Z)dZ= 



1 



dZ 



V27Tf Z—Zq 



exp- 



[ 



[log(Z-Z )- M ] : 



2f 2 



] 



(4.1.8) 
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Figure 18. Distribution of scatter signal phase angle <f> for 
various cross correlations p. 



The corresponding short-term amplitude fluctuation 
must now be considered as a conditional probability 
density; that is, given a or Z, the distribution of R is, 



P(E\zy. 



(i in; 



log 4 exp 



_{££»}. 



(4.1.9) 



the diffracted fields 50 km beyond the horizon, so 
that E must be omitted as the source of the scatter 
signals. This leaves the integral term in eq (3.2.4), 
which we have identified with single scattering of 
direct rays by the refractive blobs, as in figure 19a. 
The initial field E is that which would be received 
at the scattering blob by line-of-sight propagation 
from the transmitter. The scattered wave proceeds 
from irregularity to receiver along the second 
(dashed) line-of-sight path, and is thus generically 
similar to the relay link problem characterized by 
figure 3a. The crucial difference lies in the scatter- 
ing region V, which is now the common (wedge- 
shaped) volume defined by the intersection of tan- 
gent planes at the transmitter and receiver (see 
fig. 19a). 





McCrossen [54] has evaluated the probability that R 
exceeds a fixed level r at any instant, by integrating 
this expression over all possible values of Z, as well 
as Ry-r . 

4.2. Electromagnetic Scattering 

The basic field eq (3.2.3) for propagation through 
a turbulent medium also generates the solutions 
which describe propagation beyond the optical hori- 
zon. The single scattering (Born) approximation 
(3.2.4) is consistently used to describe scatter propa- 
gation [7, 45, 46]. E is the field which would be 

-> 
received at R if there were no refractive irregulari- 
ties, and is therefore the spherical earth diffraction 
field at points beyond the horizon. It has been 
established experimentally [47] that the scatter field 
strengths are many orders of magnitude greater than 
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Primary Line -of -Sight Field 

Scattered Field Ray Path 

- Diffracted Field Ray Path 

Figure 19. Ray paths for scatter propagation: (a) Common 
(wedge) volume for line-of-sight ray paths, and (b) scattering 
°f diffracted wave in shadow zone. 



The scattered field represented by the integral 
term in the Born approximation (3.2.4) gives rise to 
orl hogonal signal components x(t) and y(t) which are 
in- and out-of-phase respectively relative to the 
comparatively weak phase reference E . 

x+iy=:FE s =Fk/\ d 3 rG(R~r)Ae(r,t)E (r) sin X, 



F is the receiver gain factor, 



(4.2.1) 



F 2 =A R - 



* 2 G R 



4ir 



which transforms field strength (volts/cm) into re- 
ceiver voltage (volts). We omit the polarization 
factor sin x for the time being. Separating (4.2.1) 
into real and imaginary components 

x(t)=FJ <PrAe(7,t) Re[-kjG(R^)E (r)\ 

V(t)= F j v d z rAe(r,t) Im [~k}G(Rs)E (rfl 

(4.2.2) 



completes the identification of voltage components 
in figure 16a and eq (4.1.1). 

The variance and cross correlation of the orthog- 
onal voltage components required in (4.1.2) can be 
established from these expressions using the spectrum 
method, as given by (2.3.1). Interchanging the 
orders of wavenumber (k) and scattering volume (r) 
integrations yields compact results similar to (3.4.1) 
in which the turbulence model S(k) and geometrical 
considerations are explicitly separated. 

^ = ^lC (Pks ^\ f dW k ?Im[k}G(Rf)E (r)l 
(xy)= : £-Ad'kS(k)I f tPi^Be [kjG Ar)# (r)l\ 
x/ f dh^Im[kjG(R^)E (r)Tl (4.2.3) 

The last integral is real in virtue of the Hermitian 
property of the spectrum S*(k) = S(—k), which 
follows from eq (2.3.1). The two basic scatter 
propagation integrals contained in (4.2.3) have not 
been evaluated explicitly to date. They have been 
attacked by the author at two levels of sophistication. 
The first method employs the simple common volume 
cutoff indicated by figure 19a. E Q (r) is taken as a 
spherical wave expanding from the transmitter, and 
G(R,r) is the free space Green's function of eq 
(3.2.5), representing unshielded expansion of the 
scattered wave from the blob to receiver. The 
common volume is then the wedge-shaped region 
described above and it is convenient to use rectangu- 
lar coordinates centered on the wedge apex and the 
great circle plane between transmitter and receiver. 
The second approach is indicated by figure 19b 
and recognizes the role of both line-of-sight and dif- 
fracted primary and/or scattered waves. In this 
case, one takes E to be the actual (series) solution 
for a dipole radiating above a spherical conducting 
earth. This automatically includes line-of-sight 
propagation to blobs in the common volume, as well 
as blobs below the horizon which are reached only 
by the diffracted component of E . Subsequent 
propagation of the scattered radiation must also 
recognize the earth-screening effect, which is to say 
that the Green's function G(R,r) must be that appro- 
priate to the exterior of a conducting sphere. How- 
ever, this is just the field due to a unit dipole placed at 
the scattering point, and is therefore given again by a 
series solution of the form used for E . This ap- 
proach treats the three types of scatter paths shown 
in figure 19b exactly, and includes the scattering of 
diffracted waves and diffraction of scattered dif- 
fracted waves in a natural fashion. This refinement 
is probably of academic interest only, since Arons [55] 
has concluded that such effects are quite small for 
normal scatter circuits. This approach does have 



the advantage of treating height-gain effects and 
near horizon paths in a unified manner. The volume 
integrations in (4.2.3) required for this approach are 
not trivial and the surviving summations must be 
performed by the Watson transformation [56] in 
the end. 

All of the precise development described above 
is essentially new and unpublished. The con- 
siderable theoretical literature on scatter propaga- 
tion is based on the cross section approach which is 
built on a series of approximations to the integral 
representation (4.2.1) for the (complex) scattered 
field, and treats only the average scattered power 
(R 2 )=(x 2 )+(y 2 ). The first step is the far field 
approximation, which notes that the scattering 
blobs are always many, many wavelengths from 
both the transmitter and receiver (see, however, 

ref. 24). This allows the scalar distance between 

-> -> 

the receiver R and variable scattering point r to be 
expanded as follows : 



~R 2 ~Kf==R 2 -*f> (4.2.4) 



B 2 



where k 2 is the propagation vector of the scattered 

-> 
wave (\k 2 \=k f =2Tl\). The field E incident on the 
point r is representable as, 



E °-\-^r r; 6 ' 



where k x is the upgoing (line-of-sight) propagation 

vector. G T is the gain of the transmitter along 

-> 

ki, and P T is the transmitted power. Combining 

these expansions with (4.2.1) yields: 



l 



^W,^ 



R\R 2 



47T 



The mean square received power is given in terms 
of receiver aperture A R = G r \ 2 /4t by 



■\GrGrGtGt 



RIR2R1R2 



(4.2.6) 



where k\, k 2 and the gain functions are implicit j 
functions of position. The second approximation 

is now invoked by noting that the space correlation 

-> -> 

in (4.2.6) vanishes unless the points r and r r are 
within several hundred meters of one another. 
The other factors in the integrand are relatively 
insensitive functions of position, so that 

Ri=Ri, H 2 =R 2 , k\ — k 2 =ki — k 2 , G T =G T , G R =G R 
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over the important joint regions of integration, 
transforming to sum and difference coordinates, 

r=r l +r 2 and p=r x — r 2 , 



By 



one can rewrite (4.2.6) as 



, _p JL r 

*~16A 2 Jt 



dPr 



G R G 7 
R\R\ 



f (PpeW^W <Ae 2 )C(p). 



Because C(p) is nonvanishing only over a small 
volume of p space, the p integration can be extended 
over all space without misrepresentation. The 
Fourier inverse of (2.3.1) then indicates that this 
(difference) integral is nothing more than the 
spectrum S(k) evaluated at the uniquely important 
"scattering wavenumber", 



K- 



47T 



A 



(4.2.7) 



where 6 is the scattering angle formed by line-of- 
sight rays to the volume element d?r, 



"16X 2 



j; 



s 



(Pr- 



[(?**)] 



mm 



G T G R - (4.2.8) 



This expression indicates that the scattering process 
selects those blobs in the hierarchy described by S(k) 
whose wave numbers satisfy (4.2.7). This corre- 
sponds to constructive interference of the waves 
diffracted by the appropriate grating (i.e., blob) 
spacing for the Bragg angle 0. This selective isola- 
tion of only a certain size class of blobs is equivalent 
to the narrow band filter shown in figure 20 applied 
to the spectrum of irregularities. 



S(k) 




Figure 20. Typical spectrum of refractive irregularities show- 
ing narrow band filter at the scattering wavenumber which 
characterizes tropospheric scattering. 

is the average scattering angle and X is the radiation wave length. 



The expression (4.2.8) for received power is often 
quoted in terms of a scattering cross section 

<r(0,X)=£s(ysinf). (4.2.9) 

This cross section measures the energy scattered per 
unit volume V, into a unit solid angle dec which makes 
the angle 6 with the initial Poyting vector S . This 
is illustrated in figure 21. The scattering cross 
section is always defined for a unit incident power, 
so that at a distance R from the scattering element V, 



IP 



E S (R) 



E 



(4.2.10) 



which is equivalent to (4.2.9) in virtue of the Born 
approximation (4.2.1) for E s . a has the units cm -1 , 
since it is an ordinary cross section (cm 2 ) per unit 
volume (cm -3 ). 




Figure 21. Geometrical description of scattering cross section 
a per unit volume and solid angle dil per unit incident power. 



The angular dependence of a is often interpreted 
as the scattering polar diagram of an average blob. 
Reference to table 1 shows that almost all of the 
scattered radiation is projected into a small cone of 
beamwidth X// about the forward direction. This 
means that scatter propagation circuits depend on the 
extremely weak side lobes of the scattering blobs, 
which explains, in part, why large dishes and high 
power are required to operate such links. Fesh- 
bach [57] has suggested that many small angle 
multiple scatterings could be more efficient than a 
single large angle scattering. Such effects certainly 
become increasingly important at frequencies above 
1,000 Mc and/or very long paths, since the con- 
vergence of the Bom series is related to the smallness 
of the equivalent line-of-sight rms phase variation 
computed along the equivalent (kinked) ray path 
eq (3.2.11) with Lead. If this view is correct, the 
entire subject, which is summarized in subsequent 
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sections, requires revision. Further research should 
be directed to this interesting point. 

The equivalence between a and S given by eq 
(4.2.9) can be used to produce the familiar expression 
for received power P R ^(R 2 ) to transmitted power 

PtH^oI 2 



Pr 



16tt 



[<Pr 



cr(0,X) 

mm 



G T G R . (4.2.11) 



The corresponding spectrum integral is more funda- 
mental, since the joint frequency-distance dependence 
is emphasized thereby. 



Pr = J_ 
P T 16X 2 



f ^3 GtGr o /4tt . 0\ 



(4.2.12) 



Data taken with very narrow beams (see fig. 16b) is 
usually expressed in terms of this ratio. 6 If the 
beams are so narrow that the scattering angle, and 
hence the spectrum, is substantially constant over 
the common volume, the integration can be collapsed 
to 



P^ 

Pt~ 



VG T G R 



S 



/4tt . r 0\ 
U Sm 2> 



(4.2.13) 



since R 1 ^R 2 ^d/2. The scattering volume is set by 
the aerial beamwidths £ and the transmission distance 
d by F^2(/3d/4) 3 /0, where the average scattering 
angle is d=d/a for d<^.a. 
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Pt~ 



(?G T G R /2icd\ 
' 32X 2 d 2 \\aj 



(4.2.14) 



The majority of broad-beam data is given with 
respect to the free-space power P FS which would be 
received over a line-of-sight path of the same 
length d. 



G T \ 2 G E 



4wd 2 47r 



(4.2.15) 



This measure of transmission efficiency emphasizes 
the loss which is due to the scattering propagation 
agency alone, 

The receiver and transmitter gains also cancel in this 
expression, since they are relatively constant over 
the scattering volume. The equivalent expressions 
(4.2.12) and (4.2.16) will be used interchangeably 
in what follows. 



6 Additional factors of two are sometimes included [49] in such expressions to 
represent the addition of ground reflected waves. Our gain factors are defined to 
include all such effects. 



4.3. Radiofrequency Dependence 

The variation of received power with carrier fre- 
quency can be established without evaluating the 
difficult volume integrals presented in the last section. 
This is because the value of X and 6 relevant to most 
scatter paths combine to emphasize "scattering- 
wave numbers' ' which lie in the inertial range of the 
spectrum of tropospheric irregularities. 7 The physi- 
cal models for S(k) are proportional to simple inverse 
powers of k in the inertial range, and the radio- 
frequency dependence of the power ratio expressions 
(4.2.12) or (4.2.16) can be brought outside the 
integrals. If S(k)=S k~ n , one has from (4.2.12), 



\ n ~ 



16 



x 



d s r 
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Sq 
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(4.3.1) 



There is still a frequency dependence concealed in the 
gain functions, which is due to both wave length- 
dependent height gain factors and the absorbing 
area-gain relationship G R =A R 47r/X 2 . The ratio of 
received powers for scaled aerials on the same path 
is therefore 



p r M 
Pr^Y 



<® 



(4.3.2) 



The mixing-in-gradient model (2.4.6) corresponds 
to n= 5 and predicts a linear-wavelength dependence 
of scattered power. This prediction is in good agree- 
ment with careful analysis of the data from nu- 
merous broad-beam experiments when the role of 
height gain factors is carefully included [58]. The 
same model implies a X 3 variation for the ratio 
Pr/Pt, which is quite consistent with narrow beam 
data [59]. It would appear that the X~ 033 dependence 
predicted by the Obukhov model (2.4.4) is not con- 
sistent with present radio data. 

Bean [60] has studied the effect of tropospheric 
absorption on scatter signals as a function of fre- 
quency and distance, and finds that such effects do 
not influence the above conclusions. The variability 
of the frequency scaling exponent has been noted 
on all links, and one must be quite explicit about 
just which quantities are being compared on the 
two frequencies. It would seem that ratios of hourly 
means on the same path for the same hours would 
be the best choice of data, since [1] two widely 
separated frequencies do not fade together (see 
section 4.9), and [2] diurnal and seasonal effects 
may affect two frequencies differently. On the other 
hand, the NBS group finds that hourly medians on 
adjacent identical paths are not always the same, 
and a further examination of the problem is indi- 
cated. The above statements about frequency de- 
pendence refer to the average power, and hence to 
the average frequency exponent. It should be em- 
phasized that it is this average value, and not the 
extreme values, [58, 61] which should be correlated 
with present (smoothed) experimental data. 

7 For example, if /= 1,000 Mc and rf=400 miles, X=l ft and 0=0.1 radian eq (4.2.7) 
gives K=2 m-i. 
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4.4. Distance Dependence 

A large amount of experimental data on signal- 
power attenuation relative to free space drawn from 
reference [62] is plotted in figure 22. This is a hetero- 
geneous sample of long-term medians, obtained in 
most cases with broad beam antennas. The data 
were taken during the wintertime, mostly in the after- 
noon, when tropospheric scatter is expected to be 
the dominant propagation mode. These data have 
not been normalized for the effects of radiofre- 
quency, antenna heights, terrain or meteorological 
influences, which explains the considerable data scat- 
ter. If the observed linear wave-length dependence 
of the scatter signals is used to reduce these data to 
a common radiofrequency base, the data scatter is 
reduced. The empirical curve drawn through the 
data corresponds to an exponential variation of 
signal power at large distances, with approximately 
0.1 db/mile. At shorter distances, the distance 
dependence is more complicated. 

The distance dependence of scattered power is 
not as simple to extract as the radiofrequency vari- 
ation. This is because such predictions depend 
upon careful integration of the spectrum's angular 
dependence over the common volume V, which, in 
turn, is defined by the transmission distance d. The 
gain factors are essentially constant in (4.2.12), so 
that the distance dependence of the mixing-in-gradi- 
ent model (2.2.4) can be expressed in terms of the 
integral 



V^ I (Pr 



1 
R\R 



\dh) 
(sin£f 



(4.4.1) 



where all factors in the integrand depend on the great 
circle path distance d. The mixing-in-gradient model 
is subject to a more severe test than other models 
because it contains no (adjustable) turbulence pa- 
rameters. The NBS group has used the meteoro- 
logically observed exponential variation of </e () /dh 
up to stratospheric heights to predict the experi- 
mental data presented in figure 22 quite well out to 
700 miles [61, 62]. 

A substantial amount of early theoretical predic- 
tion [58, 63, 64] was based on height dependences 
assumed for the mean square dielectric fluctuation 
<Ae 2 >oc/r w . This also increased the distance scal- 
ing law exponent by n and supplied some of the 
missing exponents for older correlation function 
models thereby. Insofar as one believes in the 
physical theories of mixing models, it is more realistic 
to insert measured profiles for e (h) and let the 
theories be tested thereby. 

In making precise comparisons between theory and 
experiment, it is necessary to recognize several 
important details. The effect of (irregular) terrain 
obstacles in modifying the takeoff angles, and hence 
the minimum scattering angle versus distance rela- 
tion, led the NBS group to the concept of effective 
angular distance [58]. Refraction of the upgoing 
and downgoing beams by the mean profile e has 
been analyzed thus far by the equivalent earth radius 
method [58]. Bean [60] established theoretically the 
corrections due to tropospheric absorption on dis- 
tance dependent scaling laws. 

The simplest case to analyze is that for narrow- 
ant enna beams, since the volume integration is then 
collapsed and approximated by eq (4.2.14). The 
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22. Experimental data on distance dependence of tropospheric scatter field strength 
relative to free space for all frequencies, after Rice, Longley, and Norton [62]. 
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mixing-in-gradient model (2.4.6) predicts the follow- 
ing distance scaling law, 



p.* 
p 



t [_dh J v 



constant 



(4.4.1) 



for identical antennas and height gain effects. This 
result is independent of turbulent parameters, 
whereas the Obukhov model (2.4.4) gives 

Both models depend upon the gradient of the mean 
profile € at the average scattering altitude. This 
point rises as the path length is increased and the 
gradient probably decreases accordingly. The nar- 
row-beam data [65] suggests a scaling law of the 
form 1/d 9 , or possibly even e~* d . Appropriate 
fitting of de /dh and/or kp in the foregoing can bring 
either theoretical result into agreement with almost 
any experimental scaling law, so that this type of 
measurement cannot be considered a very severe 
test of the theory at the present time. 

4.5. Absolute Power Estimates 

Probably the most difficult problem in scatter 
theory is to make reliable absolute power estimates. 
All of the explicit and implicit factors in (4.2.12) 
must be carefully evaluated. The following points 
must be recognized in reaching accurate power 
predictions: (a) Horizon limitations of local terrain, 
(b) height gain effects in the antennas, (c) lobe 
patterns of the antennas, {d) variation of dielectric 
irregularity intensity with height, (e) refraction of 
primary and scattered rays, (f) absorption, (g) 
anisotropy influences, and (h) transmitter and 
receiver antenna loading. 

The problem has only been solved completely by 
the NBS group, using numerical techniques [61, 62]. 
They considered broad-antenna patterns with the 
mixing-in-gradient model and obtained good (ab- 
solute) agreement with experimental data. 

The problem is primarily one of finding a judicious 
coordinate system in which to perform the com- 
plicated volume integration of (4.2.12). The NBS 
spherical system [66] can be integrated exactly for 
many turbulence models [67], but is not adaptable 
to height-dependent scattering cross sections. The 
rectangular system used by Staras [64] circumvents 
this problem, but relies on approximations which 
are difficult to assess. It is anticipated that the 
precise integration schemes suggested by eq (4.2.3) 
will supply an accurate means for making absolute 
power level comparisons. A height-dependent 
factor (de /dh) 2 in the spectrum enters into (4.2.3) 
as a factor de /dh in each fundamental integral. 

An important achievement for scatter theory is 
the reliable prediction of signal strengths beyond 
the horizon from surface measurements of (de Q /dh) 
at the earth's surface [62]. This means that diurnal 
and seasonal (relative) variations of scattered power 



can be anticipated without recourse to complicated 
volume integrations, since the gradient profile, 



6?€ _r(i€o"| 
dh L^Jc 



F(h), 



(4.5.1) 



permits the time-variable ground level factor to be 
removed from the geometrical integrals. 

4.6. Amplitude Stability of Scatter Signals 

The characteristic random fading of the amplitude 
of scatter signals is of considerable practical im- 
portance, since useful information can only be 
transmitted when the signal strength is above a 
fixed level. The temporal variations of the signal 
envelope are physically related to the incoherent 
addition of (random) Doppler-shifted waves received 
from different elements of the scattering volume. 
A typical time record of scatter signals is shown in 
figure 23, together with the time-displaced vector 
voltage diagrams corresponding to the instants t\ and 
t 2 . The basic theoretical problem is to predict the 
joint probability that the signal amplitude will 
assume the values R x and R 2 at times separated 
by an interval T~t 2 —t\. The average time between 
maxima and zero crossings, amplitude auto correla- 
tions, etc., can all be computed from such a result. 

The purely statistical descriptions of fading must 
be discussed first. The distribution of time-dis- 
placed signal amplitudes shown in figure 23 depends 
on the four dimensional probability density for the 
set (x u y h x 2 , y 2 ). Since the individual orthogonal 
signal components x(t) and y(t) are considered to be 
distributed in a Gaussian manner at each instant, 
the required density is a four dimensional Gaussian 
multivariate form, 

P{Zi)= aJ\M\ eXp "~4 ? ? Z^MTj 1 (4.6.1) 




<R> 





Figure 23. Typical time record of randomly fading scatter 
signal amplitude with corresponding time-displaced vector 
voltage diagrams. 
(R) is the average amplitude and <f> is the instantaneous phase of the signal. 
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whore Zi==x(f), Z 2 =y(t), Z z =x(t+r) J Z A =y(t+T). 
of the fading signal. 

-<z 2 (0> (x(t)y(t)) 

(*(t)y(t)) (v\t)) 

(x(t+r)x(t)) (x(t+r)y(t)) 



The moment matrix M tj is the basic characterization 8 



M-- 



(x(t)x(t+r)) 
(x(t+r)y(f)) 

WV + t)) 



(x(t)y(t+r)) 

(y(t)y(t+r)) 

(X(t+r)y(t+r)) 



l(x(f)y(t+r)) (y(t)y(t+r)) W+rMt+r)) (y 2 (t+r)) 



4.6.2 



The matrix M tJ is symmetrical because the compo- 
nent processes are assumed to be stationary. The 
cylindrical transformations x=R cos <j> and y=R 
sin (l> which define the signal amplitude and phases 
(see fig. 23) must be inserted into expression (4.6.1) 
to establish the time-displaced phase and amplitude 
distribution. This portion of the problem is ana- 
lytically very complicated unless the moment matrix 
Mtj assumes especially simple forms, and it is desir- 
able to investigate the time correlations themselves 
before pursuing the statistical problem further. 

The propagation problem is to compute the 
moment matrix of eq (4.6.2) explicitly. This has 
not been done. One can, however, establish general 
expressions for the time-displaced correlations 9 in 
Mij by using the integral representations for x(t) 
and y(t) established previously from the Born 
approximation (4.2.2), viz, 



*(*)=[ <PrAe(r,t)H(r), 

y(t)= f d z rAe(r,t)J(r), (4.6.3) 



where H(r) and J(r) are the real and imaginary 
parts respectively of the complex product [k 2 G(R,r) 
E (r)] in eq (4.2.2). The appropriate averages can 
now be computed using the integral representation 
(2.5.2) for the space- and time-displaced correlation 
of dielectric fluctuations. 



(x(t)x(t + r))= f d?rH{r) f dVff(rO<A € (r,*)Ae(r',*+T)> 

J V J V 

=7^(d*kS(k)rilk,T]e-^\ f fflrH(r)fi" 

\^J J \Jv 

(x(t)y(t+r))=(x(t+r)y(t)) 

=7^( d s k S(k)rj(k,t)e-^ T (I dhH{r)^\([ (Pr f J(r , )^' r, \ (4.6.4) 



The important point to note here is that the geo- 
metrical integrals involving H(r) and J(r) are pre- 
cisely those required for the same-time calculations 
of signal properties. Although these integrals have 
not been evaluated to date, it is significant that the 
two-time problem is no more difficult than the same- 
time problem as far as the propagation calculations 
are concerned. 

The statistical problem has been treated hereto- 
fore [51, 69] by assuming very special forms for the 
moment matrix. 



(xl)=(xi)=(y{)=(yl)=°* 

(x 1 x 2 ) = (y 1 y 2 ) = a 2 f x(T) 



(4.6.5) 



/x(r) is here called the basic propagation correlation 
coefficient. All cross correlations are assumed to 



8 \M\ is the determinant of Ma and Af»yi its inverse. 

9 The one-time components are special cases of the time-displaced correlations, 



vanish identically. The joint probability density 
(4.6.1) now reduces to a comparatively simple form, 
which is readily transformed into the polar amplitude- 
phase coordinates of figure 23, 



P[R u R 2 ,(t>i,(t>i',r] 



R 1 R 2 



"47T 2 0- 4 (1- M 2 ) 



exp 



-[ 



R 2 1 +Bl-2^R 1 R 2 cos((t> l - <fe) 
2<r 2 (l- M 2 ) 



] 



(4.6.6) 



The joint amplitude distribution for B 3 and R 2 is 
obtained by integrating this expression over both 
phase angles. 
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The amplitude autocorrelation can be expressed in 
terms of complete elliptic functions, 

(R(t)R(t+r))=(R 1 R 2 )=a 2 {2E( f ,)-(l-^)K^)} 

(4.6.8) 

and is plotted in figure 24 as a function of /x(r). 
For large time separations (i.e., /x small), this result 
approaches the square of the average amplitude, 

(R) 2 =- a 2 . For small separations or /x near unity, it 

z 

approaches (7? 2 )=2cr 2 . The frequency spectrum of 
amplitude variations can be obtained from (4.6.8) 
by the Fourier cosine transformation, if /x(r) is 
known explicitly as a result of propagation calcula- 
tions. 
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Figure 24. Normalized autocorrelation of scatter signal 
amplitude versus basic propagation correlation coefficient 
li (r). 

The envelope fading rate, or average number of 
positive crossings of the median signal level R per 
unit time in figure 23, is a convenient measure of 
the time-var}dng random signal. This fading rate 
N can De identified in the small argument expansion 
of the amplitude autocorrelation [30]. 

(R(t)R(t+r))=(R*)[l-y 2 (irNTy . . .). (4.6.9) 

If the basic propagation correlation coefficient /x(t) 
has a similar expansion, 



m(t)=i-HM 2 



(4.6.10) 



Reduction of expression (4.6.8) leads to the following 
relationship : 



N= 



V2 



(4.6.11) 



The basic propagation problem is to predict the 
radiofrequency and distance dependence of 7. 

Early studies of fading rates on scatter links 
recognized that the primary fading mechanism was 
due to Doppler-shif ting of the signal by motion of the 
scattering medium. An rms sum of the drift velocity 
U and average eddy velocity V was imagined to 
give the effective speed [32, 68] 

N= f - [(217 Sin |J +| V{J , (4.6.12) 



where d is the average scattering angle in the com- 
mon volume. The linear frequency dependence is 
characteristic of simple Doppler-shifting. 

The geometry of a broad-beam scatter link shown 
in figure 25 indicates, however, that the situation is 
rather more complicated [69]. Consider first the 
midpoint (A) in the scattering volume. Since 
Doppler shifts are only caused by net elongation (or 
contraction) of the propagation path, the predomi- 
nantly horizontal drift will cause no fading of the 
signal scattered from midpoints. On the other hand, 
the smaller random velocity V has an arbitrary 
direction and can impose Doppler shifts on this 
component of the scattered signal. Consider next 
an off center scattering element (B). The drift 
velocity U now has a small component U± which 
bisects the propagation paths, and therefore does 
impose a Doppler shift. The fading contribution of 
these off-axis blobs is small because: (i) U± is a 
small projection of U, and (ii) the power scattered 
from point B is considerably smaller than that from 
point A due to the larger scattering angle 0'. Direct 
confirmation of this prediction was obtained by the 
Lincoln Laboratory Group [59], who noted that the 
fading rate decreased when very narrow-antenna 
beams are swung simultaneously away from the 
path midpoint A. 

A considerable advance in analyzing this fading 
phenomenon was made recently by using the spectral 
representation of the turbulent dielectirc fluctuations 
[14]. The explicit separation of drift and random 
velocity contribution to the space-time correlation 
of Aein eq (2.5.2) suggests the following time-displaced 
generalization of the cross section result (4.2.12), 



Trans. 




Figure 25. Geometrical description of doppler shifts imposed 
-> -> 

by horizontal drift speed U and random velocity V at:- (A) 
Path midpoint, and (B) off axis scattering volume elements. 
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(R(t)R(t+r)c 



A <PrS(\K\) e^'i, [Vok 1/3 \K\ 



2/3 „ 



(4.6.13) 

where K is the scattering difference vector defined by 
(4.2.7). Since the major component of the drift 
speed is usually horizontal, and therefore perpendicular 

to K for midpoint scattering, the significant contri- 
bution of drift is made by volume elements far from 
the midpoint. The random self motion V has no 
such projection factor and contributes throughout 
the scattering volume. Expansion of (4.6.13) in a 
power series of r and identification with expression 
(4.6.9) indicates that the fading rate due to self- 
motion is proportional to / 2/3 , while drift motion 
gives the usual linear Doppler dependence on fre- 
quency. Available data indicates a variable fre- 
quency exponent less than one [32], so that some 
mixture of the two effects is indicated. Further 
experimental results are needed, especially data on 
the distance-dependence of N. The theory needs 
to be clarified by evaluating n(r) from eq (4.6.4) for 
several propagation models. 

4.7. Phase Stability of Scatter Signals 

The phase stability of scatter signals has received 
little attention to date — either theoretically or 
experimentally. The question is potentially quite 
interesting because of the possibility of passing stable 
frequency references beyond line-of-sight by this 
mode. Frequency excursions about a single carrier 
are completely equivalent to time rates of change of 
the signal phase <j>(t) defined in figure 23. The mean 
frequency change is thus given in terms of the auto- 
correlation of phase 



(8f 2 ) = (4>(t)<i>(t + r)) 



= -£<*(*)*(*+t)>| 



(4.7.1) 



The relative change in phase between two instants, 

(!0(/)-0(^ + r)| 2 )-2[(0 2 )-(0(O«a + r))], (4.7.2) 

is independent of the absolute path phase reference, 
and is probably easier to measure experimentally. 

The central quantity in such studies is the auto- 
correlation of phase taken between two instants. 
For the special assumptions indicated by (4.6.5), it is 
possible to derive this result in terms of the basic 
propagation correlation coefficient jli(t). Integrating 
(4.6.6) over the amplitudes R x and R 2 gives the dis- 
tribution for </>! and </> 2 deduced by Bunimovich [70] 



P[<t>l,<t>2,T] 



,(1-M 2 ) 



4tt 2 



n i • -i 
^+ sin l co 

(1-co 2 ) 3 ' 2 . 



(4.7.3) 



where co=/x cos (</>i — </> 2 ). Wagner 10 has used this 
result to calculate the phase autocorrelation, 



'R.J. Wagner, to be published. 



(4>(t)4>(t + r))=^sm-^ 

+}- [sin "VP-tt [ dx[an- l (n cos .r)] 2 (4.7.4) 

" ,/ 

and this is plotted in figure 26 versus /-i(r). For 
very large r or /x small, the phases are completely 
uncorrelated. For /x unity, expression (4.7.4) ap- 
proaches the mean square value (</> 2 )=7r 2 /3. 

It is interesting that the phase correlation can be 
inferred from experimental data on amplitude corre- 
lations, since /x is only a parameter which may be 
eliminated between eq (4.6.7) and (4.7.4). When 
the more general moment matrix (4.6.2) is used, 
however both results will depend on two or more 
basic functions of t. In that case, it is still possible 
to convert phase into amplitude data and vice versa 
by eliminating r, although one must first complete 
explicit progagation calculations for mM> etc. 

4.8. Space Correlation and Diversity Reception 

A combination of two or more spaced receivers is 
frequently exploited to overcome signal loss during 
deep fades. This diversity reception depends upon 
spacing the antennas sufficiently far apart to insure 
independent fading of the two signals. The best 
experimental data available is that taken at 1,046 
Mc by the NBS group [72]. The statistical problem 
is very similar to that illustrated in figure 26, if 
one reinterprets Ri as the instantaneous signal 
amplitude induced in the first receiver and R 2 that 
induced in the second, at the same instant. The 
joint probability density for R x and R 2 is given by 
eq (4.6.7), if the basic correlation coefficient u(D) 
is now the space correlation between the orthogonal 
signal components at the two receivers. Staras [73] 
computed the probability that the greater of the two 
(diversity) signals is above a specified level, as a 
function of the correlation coefficient /x (D). These 
statistical predictions have been checked against 
experimental data [74] and seem to fit quite well. 

The propagation problem is to estimate the space 
coefficient u (D). Gordon [63] and Rice [68] pre- 
sented qualitative analysis which showed that the 
correlation distance normal to the propagation path 
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Figure 26. Autocorrelation of scatter signal phase versus basic 
propagation correlation coefficient /jl (r). 
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(horizontal and vertical) ought to be of the order 
\a/d, where d is the transmission distance and a 
the radius of the earth. A correlation distance 
\a 2 /d 2 along the path was also estimated. The 
RCA experiments [74] at 390 Mc for path lengths 
of 73 and 153 miles, however, failed to show the 
predicted decrease of horizontal and vertical corre- 
lation lengths with path distance. On the other 
hand, the simple wavelength dependence of these 
results is borne out by multiple frequency experi- 
ments [53, 75]. 

Analytical propagation calculations of all three 
correlation functions 11 as functions of receiver 
separation distance were derived by Staras [64], 
using the exponential space correlation model. 
His results have the same form as Gordon's, but are 
numerically different by substantial factors. Staras 
also included the effects of anisotropy in the refrac- 
tive irregularities and found that this affected the 
predicted horizontal (transverse) correlation distance, 



S-- 



1.4ir/aV 



(4.8.1) 



where r is the vertical-to-horizontal scale length 
ratio. By fitting data taken on 100 and 1000 Mc 
[53, 72] he suggested that this ratio should be about 
0.25. Such comparisons are important in that they 
bear directly on the question of anisotropy, but more 
good data is needed. 

The realization of plane-earth height-gain advant- 
ages is limited by space coherence between the actual 
receiver and its image. The usual height-gain curve 
is realized if the antenna is less than one-half the 
vertical correlation distance (\a/2d). Only slight 
improvement is gained above this height by the 
gradual decrease of the scattering angle. 

4.9. Frequency Correlation and Medium Bandwidth 

The signal-frequency bandwidth which the scatter 
medium can support without serious distortion is 
intimately related to the multipath delays which 
are experienced. These depend upon the size of 
the effective scatter volume, which in turn is de- 
termined by the antenna bandwidths and/or the 
scattering pattern of the blobs. When the antenna 
patterns (jS) are broad, the multipath is controlled 
by the scattering blobs and the approximate result 
of Gordon [63] predicts that 



A/= 



5ca 2 



*»i 



(4.9.1) 



or numerically 

A/ Mc =30/d 3 , 
where d is in hundreds of miles. This result shows 



ii Along, horizontally normal (transverse), and vertically normal to the 
propagation path. 



that the bandwidth decreases rapidly with distance. 
Rice's estimate [68] of A/ varies inversely with the 
path length d, but is based on an oversimplified 
physical model. 

A real improvement in bandwidth capabilities is 
obtained with very narrow beams. Although some 
antenna gain degradation is experienced with 
narrow beams, the angular distribution of multipath 
signals which are admitted to the receiver is reduced 
and the bandwidth thereby increased. Booker and 
de Bettencourt [76] examined this problem for beams 
narrower than the average scattering angle d=d/a 
and found 



A - Aca 



fi« 



d 



(4.9.2) 



This result indicates steady improvement with an- 
tenna size and decreases less rapidly with distance 
than (4.9.1). When aircraft fly through the path, 
the solid bistatic reflections from them supply a 
multipath signal of unusual magnitude and the above 
estimates for broad-antenna beams are meaningless. 
Figure 27 indicates typical successive plots of 
amplitude versus frequency, such as might be ob- 
tained by sweeping a narrow-band transmitter 
about the carrier f . The joint amplitude distribu- 
tion (4.6.7) can be used to describe this situation if 
Ri is interpreted as the signal strength at f Q and R 2 
that at/o+A/. The basic correlation function /x is 
now a function of the frequency separation A/ and 
other propagation variables. If the correlation is 
high, it means that the signal components at either 
end of the frequency band being used are behaving 
essentially in unison, and little or no distortion is 
expected. Staras [64] evaluated substantially this 
function directly for omnidirectional antennas with 
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Figure 27. Typical successive amplitude records plotted versus 
carrier frequency over a band A/. 



The receiver threshold level indicates the role of absolute power level in correla- 
tion comparisons. 
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the exponential space correlation model, and finds 
that the frequency correlation drops to 0.5 when 



&JntC8 73 



(4.9.3) 



where d is in hundreds of miles. Keference to 
figure 24 indicates that (4.9.3) refers also to the 
envelope (R) correlation, so a factor of eight differ- 
ence exists between (4.9.1) and (4.9.3). This prob- 
lem evidently requires further clarification and care- 
ful comparison with experimental data. It is 
possible to use the spectral separation of turbulence 
and propagation effects indicated by (4.6.4) to 
calculate the frequency correlation coefficient m(A/). 
One need only evaluate the results of the frequency- 
dependent geometric integrals at displaced fre- 
quencies and set r=0 in (4.6.4). 

The published experimental data has been taken 
with narrow-antenna beams, and is therefore not 
directly useful in resolving the theoretical conflict. 
The recent Lincoln Laboratory data [77] measures 
the frequency correlation in the manner of figure 
27, by using a rapid frequency-change system. 
Their results indicate that signals separated by 2 
Mc at 2,300 Mc over a 188-mile path give an ampli- 
tude correlation of 0.4. This is consistent with 
the television experiments of Tidd [78] on a similar 
path at 5,000 Mc. 

It is important to note that the receiver threshold 
level influences the apparent bandwidth of a system. 
Figure 27 indicates how amplitude fading at time 
t 2 can reduce the frequency band over which the 
signal is correlated with the proceeding signal. 
This has lead to the concept of "instantaneous 
bandwidth/' which finds application in designing 
efficient modulation systems to work with the scatter 
propagation mode. 

4.10. Gain Loss or Antenna-to-Medium Coupling 

When very narrow beams are used on scatter 
circuits, it is found that free-space antenna gains 
are not realized. This is known as gain loss or 
antenna-to-medium coupling loss and arises from 
the fact that signals arrive at the receiver from an 
extended scattering volume V. Narrowing the 
antenna beam eventually reduces this common 
volume faster than the aerial gain is increased (see 
fig. 16b) so that a relative gain loss is incurred. 

Booker and de Bettencourt [76] analyzed this 
problem for similar conical antenna patterns at 
the transmitting and receiving ends of the circuit. 
They argued that no gain loss should be incurred so 
long as the beamwidth /3 is greater than %d/a. When 
the beams are sharper than this average scattering 
angle, they suggest that a relative power loss factor 



L=0.43 



W : 



(4.10.1) 



will be incurred. They point out that this result is 
probably too large because the height dependence 



of (Ae 2 ) was neglected and sharp beams were em- 
ployed. The calculations by Nortion, Rice, and 
Vogler [58] used an exponential correlation function 
and an inverse square height dependence of (Ae 2 ). 
Staras [79] treated anisotropic irregularities and 
dissimilar antenna patterns. His results are between 
5 and 8 db less severe than these predictions by 
Booker and de Bettencourt, and seem to give rather 
good agreement with experimental data. 

Hartman and Wilkerson [80] have reinvestigated 
the problem, using an exponential decrease of iso- 
tropic turbulent fluctuations with height. Their 
calculations produce gain loss estimates which are 
several decibels less than those predicted by Staras, 
and in good agreement with considerable experi- 
mental data. 

4.11. Beam Swinging 

It is found experimentally that the scatter signal 
amplitude decreases when narrow beamwidth an- 
tennas are jointly swung off the great circle path, 
or simultaneously elevated as shown in Figure 28. 
For azimuth swinging, the effect is simply to increase 
the average scattering angle 6 between the ray paths 
(beams) [81]. This reduces the scattered power 
because of the large exponent of sin 0/2 in the various 
scattering cross section expressions (see sec. 4.2). 

Booker and de Bettencourt [76] analyzed the effect 
of swinging the transmitting and receiving aerials 
simultaneously in azimuth and elevation. They 
used a height-independent exponential correlation 
to make estimates of this power reduction, which 
are apparently in good agreement with available data. 

The dependence of power loss on elevation swing- 
ing must also recognize the observed decrease of 
turbulent intensity with height. Since (Ae 2 ) prob- 
ably decreases exponentially with height, the con- 
comitant scattering height change Ah can effect the 
power levels significantly if it is comparable with the 
atmospheric scale height #=22,000 ft. The prob- 
lem evidently needs further analysis and comparison 
of theory and experiment. 




Figure 28. Geometry for vertical bea?n swinging experiment 
showing increased scatter angle and higher common volume. 



231 



The preparation of this paper is due in large part 
to the hospitality extended by the Boulder Labora- 
tories of the NBS and to invaluable discussions 
with my colleagues there. W. Hartman and E. 
Barrows of that laboratory were particularly helpful 
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